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This paper presents the solution for the supplier’s profit maximization 
problem with unit commitment decisions participating in single side auction 
markets of a deregulated power system. The bids from market participants 
are received by a central pool mechanism and the Market Clearing Price 
(MCP) for energy and spinning reserve is fixed. The bid quantities are 
optimized using Differential Evolution (DE) algorithm. The supplier aims to 
achieve (more) profit than that of the rival’s participating in the competition. 
A GENCO with 6-unit participating in 24-hour day ahead energy and 
spinning reserve market is used to illustrate the methodology. The bidding 
parameters of rival’s participating in the competition are calculated by multi- 


variant Probability Density Functions (PDF). The results of the proposed 
(GENCOs) y y prop 

Independent System Operator methodology are compared with Refined Genetic Algorithm (RGA). 
(ISO) Numerical results illustrate the effectiveness of the method in solving the 


Market Clearing Price (MCP) supplier profit maximization problem. 
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1. INTRODUCTION 

Ancillary services are those functions performed to support the basic services of generating 
capacity, energy supply and power delivery. Ancillary services are required for the reliable operation of the 
power system [1]. The general approach for pricing ancillary services within competitive electricity markets 
is based on fixed contracts for a certain time period between the ISO and the market participants that are able 
to provide the required ancillary services. Thereby, the ancillary services are divided into different services 
such as spinning reserves, non-spinning reserves, and Automatic Generation Control (AGC), replacement 
reserves, voltage support, and black start. The first four services can be procured by the ISO by means of 
daily competitive auction, whereas the last two services are more suitable for purchases based on long-term 
contracts [2]. In deregulated markets, there are separate auctions for each category of reserves. In this paper, 
it is assumed that GENCOs participate in the energy markets as well as in the ancillary service (only spinning 
reserve) auction markets. 

The importance of spinning reserve which is required for system reliability is focused in this paper. 
Spinning reserve is the ability of an on-line generator (load) to increase (decrease) its output (consumption) 
in a short period of time. The time period will be determined by the system but for smaller systems the time 
period is generally smaller in order to avoid large frequency deviations [3]. 

Suppliers (GENCOs) and consumers (DISCOs) participate in the bidding process of double side 
auction markets in order to maximize the profit of suppliers and benefits of the consumers. This is achieved 
by differential evolution and dealt only the energy markets and not considered the reserve markets [4]. 
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A method of building an optimal bidding strategy under market price uncertainty using information gap 
decision theory (IGDT) has been presented. A single thermal unit participating in day ahead energy markets 
without considering the reserve markets are addressed [5]. The supplier profit maximization problem is 
solved as multi objective optimization problem by considering the rival bidding and profit functions also by 
using Genetic Algorithm (GA) [6]. The supplier (decision maker) optimization problem is formulated under 
single side auction energy markets (without considering the spinning reserve markets) and their bid quantities 
are optimized using Self adaptive Differential Evolution (SaDE) [7]. 

In all the above literature mentioned, considerations have been made only to set up the energy 
markets and not on reserve markets. In this paper, an approach for providing this ancillary service is to set up 
reserve markets which run sequentially following the energy and transmission congestion management 
markets is dealt with. In this context, DE is employed to choose a GENCO’s optimal bidding strategy among 
the sets of discrete bids. 

This paper is organized as follows: Section 2 presents the market clearing mechanism in competitive 
energy and reserve markets, section 3 problem statement of profit maximization of decision makers 
submitting bids to market operators, section 4 deals with the solution methodology to find the optimum 
schedules of the supplier, section 5 presents the results and discussions, and section 6 concludes. 


2. ENERGY AND SPINNING RESERVE MARKETS 

In the restructured power systems, GENCOs will submit bid curves to the ISO, and then ISO clears 
market after collecting bids. In the ISO’s market clearing model, ISO dispatches generating units in order 
from lowest to highest bid as needed to meet demand while considering network constraints. Once the energy 
market is cleared, each generating unit will be paid according to pricing mechanism of market. Generally 
there are two pricing mechanisms: pay-as-bid and uniform pricing [8]. Under the pay-as-bid pricing structure, 
every winning generating unit gets its bid price as its income. Under the uniform pricing structure, the bid 
price of the last dispatched unit sets the market clearing price, then all units dispatched receive the same 
MCP. In this paper, the uniform pricing structure is utilized. The formation and operation of energy markets 
were discussed in [7]. Spinning reserve, a generation based ancillary service can be made competitive and 
different from energy market. The spinning reserve service can be procured by ISO through daily competitive 
auctions. The i-th supplier spinning reserve bidding function can be represented as [9] 


SP (RP )J=0; +0 RP t=1, 2, ...T (1) 
NIRO SSR) (2) 
Rimin < Rit < Rimax (3) 


where SR, is the reserve at hour t, oP Pare the intercept and slope of the spinning reserve bidding curve of 


the suppliers respectively, R; is the reserve generation output within the set of reserve limits RiminandR imax- 
The MCP for spinning reserve (MCPR) is calculated as 


Peal 2eal (4) 
The spinning reserve dispatch by each supplier can be calculated as 


Ri, = (MCPR, — 9,”)/6{Pi = 1, 2...N (5) 


3. SUPPLIER PROFIT MAXIMIZATION PROBLEM FORMULATION 
The profit maximization objective of suppliers participating in energy and reserve markets and 
competing with the other suppliers can be stated as 


Maximize: PF= RV-TC. 


Profit (PF) is defined as the revenue (RV) from the sales of energy and reserve minus the total 
(production) cost (TC). 


PF = YN, YZ, [MCP,.P;, + MCPR;.Ri, — CPt Ri] Xie (6) 
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Where Xit the ON/OFF status of the suppliers decided by ISO and Ci(Pit) is the fuel cost function of 
the suppliers. 

The constraints included are a) Power balance and b) Minimum and maximum capacity 
limits of suppliers. 
a) Power balance constraints 

The total generation (including spinning reserve) of GENCOs participating in the electricity markets 
may be greater than or equal to the demand profile of the customers. 


Der (Pit Rid) Xi = Dica De t=1,2,...T (7) 
b) Minimum and maximum capacity limit constraints 


Generation units have lower and upper production limits that are directly related to the generator 
design. These bounds can be defined as a pair of inequality constraints 


Pimin S Pit = Pimax (8) 
Rimin Š Rir < Rimax (9) 
PitRit < Pimax (10) 


4. SOLUTION METHODOLOGY 
The solution methodology of the decision maker profit maximization with optimized strategy 
problem is given as follows: 

1. Initialization and creation of parent population: Set iteration count as 1. One of the bidding parameter of 
the suppliers in energy markets (£;) and reserve markets (g;) are optimized using a suitable algorithm. 
Here DE is employed. The bidding coefficient of suppliers in energy markets (a;) and reserve markets 
(pì) is kept fixed as the cost coefficients b; and 0.5); respectively. 

2. Calculation of bidding coefficients of rivals’: The bidding parameters of the rivals’ can be determined by 
statistical approach as given below. 

The bid coefficients of rivals participating in energy markets, Bi and ai (i = 1, 2... N) obey 

a multi-variate normal distribution with the PDF given in [10] and can be expressed in compressed form as 





(a) (a)y? PoP 
(a® BOJN His (a) puo (11) 
uP (a) (8) (oi? 
‘pie it Fit it 
wherep;,, is the correlation coefficient between at ? andg! () u, u€ ) o o and of are the parameters of 
the multi-variant normal distribution. 
The bid values of rival suppliers in energy markets are estimated as 
postaa = 1.2 x 2a; 
4a = 0.15b,,40\ = 0.154; (12) 


Pit = —0.1 


The rivals’ are Taar to bid 20% above operating cost. The mean and standard deviation of a;andf;are 
specified as [ue -= 40 Aho + 4a Jana [ue L> 40? i 4 40? | respectively, with the probability of 
0.999. 


Similarly the bid coefficient of the rivals participating in spinning reserve markets, g,;andd,can be 
determined as given below. 
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MA = Oe oi = 0.5u” 


i,t 


0 = 0505 = 0.50 (13) 


Lt 2 Ut 
Vit Pit 

It is assumed that the rival bidding coefficients are same for all the 24 hours. It is not the case in practical 

situations. But in real time, the subsequent hour bids are estimated using the previous hour bidding data. 

3. Calculation of MCP and MCPR: The MCP and MCPR are calculated with the bidding data of suppliers 
and rivals’. Based on the market price, Pand Rare calculated and limit values are checked. 

4. Determination of unit ON/OFF status: IfP;,+Ri<Pimin,thenX;, = OelseX; = 1. Thus the unit ON/OFF Xit 
status can be calculated by taking an account of the constraints to be satisfied in all trading periods. 

5. Economic Dispatch: With the calculated Xit, the optimal dispatch of power Pit and spinning reserve 
power Rit are calculated using Quadratic Programming (QP). The revenue generated and fuel costs spent 
are determined. 

6. Calculation of fitness: The fitness is calculated as per equation (6). 

7. Stopping criteria: The steps from 1 to 6 are repeated until the specified maximum number of iterations is 
reached. 


5. RESULTS AND DISCUSSION 

To illustrate the optimal bidding strategy, a GENCO with six suppliers are considered to be 
participating in 24 hour day-ahead electricity market. The results of test systems with and without optimized 
bidding strategies are tabulated. The generator and load data of the test system are taken from [11] and given 
in Appendix as Table Al and A2 respectively. The proposed methodology is implemented on INTEL core, 
i3 processor, 3GB RAM and simulated in MATLAB 7.10 (R2010a) environment.The working algorithm 
used here is the seventh strategy of DE i.e. DE/rand/1/bin in which DE represents differential evolution, rand 
is any randomly chosen vector for perturbations, 1 represents the number of difference vectors to be 
perturbed and bin is the binomial type of crossover used. The seventh strategy is the most successful and 
widely used in optimization problems such as emission constrained economic dispatch [12], optimal power 
flow [13] and optimal design of gas transmission network [14]. 


5.1. Parameter Selection 

The results are sensitive to algorithm parameters. Hence, it is required to perform repeated 
simulations to find the suitable values for the parameters. Optimal parameter combinations are 
experimentally determined by conducting experiments with different parameter settings. The following 
control parameters have been chosen for the test system. 

Population size NP = 250 

Crossover Ratio CR = 0.9 

Differentiation or mutation constant F = 0.5 

Maximum number of iterations, MAXITER = 200. 


5.2. Optimal Bidding Strategy Using DE 

The supplier-6 aims to maximize its own profit and other generators (1 to 5) are its rivals’. This 
example system is utilized for 24 hour (load) demand. The spinning reserve is maintained as 10% of the 
demand. The fuel cost equation is expressed in quadratic form as 


Ci (PatRi)=a; (Pt Rit)? +b; (Pat Rit) +Ci (14) 


The supplier who is aware of market power in deregulated market is likely to bid above the marginal 
production cost. Hence the optimum values of /; and @¢ are searched using DE within the intervals [1.05x2a,, 
1.35x2a,] and 0.5x[1.05x2a,, 1.35x2a,] respectively. The algorithm used here is the seventh strategy of DE 
i.e. DE/rand/1/bin. 

Based on the optimized bidding value of supplier-6 obtained from DE technique and rivals’ bidding 
value from PDF, MCP and MCPR are fixed by PX and ISO respectively. The power dispatch and spinning 
reserve allocation of the suppliers is calculated with MCP and MCPR values in all trading hours. If the 
supplier is not able to provide minimum power requirement, then the corresponding supplier is not allowed to 
participate in the competition. Thus the ON/OFF commitment of the suppliers determined by pool operators 
for all the 24 hours and the economic power dispatch in all the trading hours is calculated using QP. It is 
observed that in the first 3 hours and 5" hour, sixth supplier is in OFF condition because of bidding scheme. 
Since in these trading hours, supplier-6 cannot be able to supply even minimum requirement. So, the 





Optimal Bidding and Generation Scheduling of Utilities Participating in ... (B. Rampriya) 


50 o ISSN: 2252-8814 


supplier-6 is made OFF during these hours. Suppose if the pool operator allows supplier-6 to enter into 


competition, then economic loss may occur during these hours. 


The commitment schedule and power dispatch in energy and spinning reserve markets of the 
suppliers are presented in Table 1. The values of MCP of energy and spinning reserve, revenue generated, 


cost spent obtained for all the 24 trading hours are tabulated in Table 2. 


Table 1. Power Dispatch in Energy and Spinning Reserve Markets 








Hour Power (MW) Reserve (MW) 

Unitl Unit2 Unit3 Unit4 UnitS Unit6 Unitl Unit2 Unit3 Unit4 Unit5 Unit6 
1 69.40 54.72 52.06 72.40 75.41 0 2.82 7.38 11.02 5.86 8.87 0 
2 73.56 57.68 54.62 76.56 79.57 0 3.32 7.71 11.31 6.33 9.33 0 
3 77.72 60.64 57.19 80.72 83.73 0 3.79 8.04 11.59 6.79 9.80 0 
4 72.64 57.03 54.05 75.64 78.64 40.00 3.39 7.16 11.35 6.40 9.40 3.70 
5 78.76 61.38 57.83 81.76 84.77 0 3.90 8.12 11.66 6.90 9.91 0 
6 76.80 59.98 56.62 79.80 82.80 40.00 3.81 8.05 11.61 6.81 9.81 3.89 
7 86.58 66.95 62.65 89.58 92.59 42.66 4.85 8.80 12.25 7.86 10.86 4.38 
8 103.47 78.96 73.06 106.47 109.47 50.57 6.73 10.13 13.41 9.73 12.73 5.26 
9 125.98 94.98 86.95 128.99 131.99 61.11 9.23 11.91 14.95 12.23 15.24 6.43 
10 135.36 101.66 92.73 138.37 141.37 65.51 10.27 12.65 15.59 13.28 16.28 6.92 
11 135.36 101.66 92.73 138.37 141.37 65.51 10.27 12.65 15.59 13.28 16.28 6.92 
12 127.86 96.32 88.11 130.86 133.86 61.99 9.44 12.06 15.08 12.44 15.45 6.53 
13 124.11 93.65 85.79 127.11 130.11 60.23 9.02 11.76 14.82 12.03 15.03 6.33 
14 125.98 94.98 86.95 128.99 131.99 61.11 9.23 11.91 14.95 12.23 15.24 6.43 
15 131.61 98.99 90.42 134.61 137.62 63.75 9.86 12.35 15.34 12.86 15.86 6.73 
16 129.74 97.65 89.26 132.74 135.741 62.87 9.65 12.20 15.21 12.65 15.65 6.63 
17 122.23 92.31 84.63 125.23 128.24 59.35 8.81 11.6 14.69 11.82 14.82 6.24 
18 109.10 82.97 76.53 112.10 115.10 53.20 735 10.57 13.79 10.36 13.36 5:55 
19 112.85 85.64 78.85 115.85 118.85 54.96 7.78 10.87 14.05 10.77 13.78 5:75 
20 118.48 89.64 82.32 121.48 124.48 57.60 8.40 11.32 14.44 11.40 14.40 6.04 
21 110.97 84.30 77.69 113.97 116.98 54.08 7.56 10.72 13.92 10.57 13.57 5.65 
22 99.71 76.29 70.75 102.72 105.72 48.81 6.31 9.83 13.15 9.32 12.32 5.07 
23 84.70 65.61 61.49 87.71 90.71 41.78 4.64 8.65 12.12 7.65 10.65 4.28 
24 74.72 58.51 55.33 71.12 80.72 40.00 3.60 7.90 11.48 6.61 9.61 3.80 





Table 2. Summary of Results of Supplier-6 








Hour Pasan Br MCP ($) MCPR($) Total Revenue($) Total Cost ($) 

1 323.99 35.95 4.3207 1.0345 0 0 

2 341.99 38.00 4.4596 1.0414 0 0 

3 360.00 40.01 4.5984 1.0483 0 0 

4 378.00 42.00 4.7373 1.0551 193.39 191.36 
5 364.5 40.49 4.6332 1.0500 0 0 

6 396.00 43.98 4.8761 1.0620 199.18 192.19 
7 441.01 49.00 5.2233 1.0791 227.54 205.88 
8 522.00 57.99 5.8481 1.1099 301.56 246.99 
9 630.00 69.99 6.6812 1.1510 415.71 308.43 
10 675.00 74.99 7.0283 1.1681 468.48 336.25 
11 675.00 74.99 7.0283 1.1681 468.48 336.25 
12 639.00 71.00 6.7506 1.1544 426.01 313.89 
13 621.00 68.99 6.6118 1.1476 405.52 303.02 
14 630.00 69.99 6.6812 1.1510 415.71 308.43 
15 657.00 73.00 6.8895 1.1613 447.00 324.96 
16 648.00 71.99 6.8200 1.1579 436.45 319.40 
17 611.99 67.98 6.5423 1.1442 395.46 297.66 
18 549.00 60.98 6.0564 1.1202 328.44 261.65 
19 567.00 63.00 6.1952 1.1270 346.97 271.67 
20 594.00 66 6.4035 1.1373 375.70 287.11 
21 557.99 61.99 6.1258 1.1236 337.653 266.63 
22 504.00 56 5.7092 1.1031 284.25 237.49 
23 432.00 47.99 5.1538 1.0757 219.93 201.57 
24 387.00 43 4.8067 1.0585 196.288 191.77 


Total Profit / day = $1487.1 





Figure 1. shows the variation of MCP and MCPR with respect to trading hours from | to 24. The 
MCP is fixed based on the load demand profile. For the given load profile, MCP is increasing from 1 to 11. 
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At 12" hour, load demand is decreased and thus there is reduction of energy and reserve price. It is observed 
that the energy price is about four times higher than reserve price. 
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Figure 1. Variations of MCP and MCPR 


The performance and effectiveness of the proposed methodology is examined in comparison to the 
solutions given by RGA [11]. Table 3 shows the comparison of hourly profit of the proposed DE with 
optimized bids and RGA methods. 


Table 3 Comparison of Hourly Profit by DE and RGA Methods 








Profit ($) Profit ($) 
Hom RGA [11] DE Hou RGA [11] DE 

1 0 0 13 93.78 102.5028 
2 0 0 14 97.94 107.2823 
3 0 0 15 110.8 122.0424 
4 2.92 2.0362 16 106.5 117.0521 
5 0 0 17 89.67 97.7936 
6 8.22 6.9913 18 62.51 66.7957 
7 22.47 21.6651 19 69.98 75.3010 
8 51.73 54.5647 20 81.62 88.5858 
9 97.94 107.2823 21 66.22 71.0132 
10 119.6 132.2337 22 44.83 46.7619 
11 119.6 132.2337 23 19.51 18.3609 

12 102.1 112.1321 24 5.54 4.5137 
, RGA 1373.48 

Total Profit ($) DE 1487.1 





It is clear that the supplier-6 receives high profit in all the trading hours when it submits the bids by 
optimizing the bid coefficients using DE to ISO. The profit distribution of supplier-6 is high in all the trading 
hours and can be seen in Figure 2. There is a net profit difference of $113.62 for supplier-6 between RGA 
and DE. 
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Figure 2. Comparison of Hourly Profit of Supplier-6 
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The convergence characteristic of the system is shown in Figure 3. The maximum number of 
iterations is fixed at 200 and the algorithm is able to converge before 150 iterations in all the runs. The best 
solutions found for the problem are tabulated. The highest profit obtained in simulation using DE is $1487.1. 
Out of 50 simulation runs, the proposed algorithm produces feasible solutions in 21 runs with best profit 
value of $1487.1. The simulation results obtained in 50 independent runs using DE are given in Table 4. 
The best and worst profits obtained are $1487.1and $1382.5 respectively. 


1500 
1480 
1460 
1440 - 
1420 
1400 4 
1380 + T T 
0 25 50 75 100 125 150 175 200 
NO. OF ITERATIONS 


PROFIT ($) 


Figure 3. Convergence Characteristics 


Table 4 Summary of Simulation Results for 50 Independent Trial Runs 








No. of trials Revenue ($) Cost ($) Profit ($) a hip at SD of profit ead 
4 8668.2 7280.6 1387.6 
a 7602.9 6128.6 1474.4 
8 8707.6 7325.1 1382.5 1442.2 52.38 1487.1 
10 7468.2 5988.9 1479.3 
21 6889.721 5402.6 1487.1 





Experimentation has also been carried out by calculating Pe and ọ6 with joint PDF and without 
optimizing these bid coefficients. In this case, ISO fixes the price and arrives at a schedule in such a way that 
the supplier-6 is put OFF during hours 1 to 5. So, the total profit of supplier-6 in this case is reduced to 
$1369.30. Table 5 shows the comparison of simulation results of optimal bidding strategy using RGA, DE 
and bidding values by joint PDF. The average execution time for a single run using DE is about 12.64 
seconds. The main difference in constructing better solutions is that RGA relies on crossover while DE relies 
on mutation operation and thus DE algorithm faces a promising approach for solving supplier optimization 
problem. 


Table 5 Comparison of Results of 6 Unit 24 Hour System 








Method Profit ($) Number of trial runs 
RGA [11] 1373.48 -- 
PDF (without bid optimization) 1369.30 -- 
DE (with bid optimization) 1487.1 21 





6. CONCLUSION 

The methodology using DE is proposed to determine optimal bidding strategy for a GENCO in 24- 
hour energy and reserve markets. GENCO submits 24 hourly supply-bidding curves for energy to the PX, 
and 24 hourly supply-bidding curves for reserve to the system operator. Based on bidder information, load 
demand and reserve, the energy and reserve awarded to each bidder are determined. The proposed method is 
developed based on the viewpoint of the GENCO as a supplier wishing to maximize profit. Investigation 
reveals that DE performs much better than GA in terms of convergence rate, quality of solution and success 
rate. The DE algorithm can solve the problem efficiently and accurately. 


IJAAS Vol. 7, No. 1, March 2018: 46-53 


IJAAS ISSN: 2252-8814 o 53 




















7. APPENDIX 
Table A1. Generator Data of Test System 
Unit ai ($/MWh’) bi ($/MWh) ci ($/h) Pmin (MW) Pmax (MW) 

1 0.0125 2 0 40 160 

2 0.0175 1.75 0 40 140 

3 0.02 1.5 0 30 120 

4 0.0125 1.9 0 40 170 

5 0.0125 1.8 0 40 180 

6 0.0275 1.85 58 40 100 

Table A2. Load Data of Test System 

Hour 1 2 3 4 5 6 7 8 9 10 1 12 
Load(MW) 360 380 400 420 405 440 490 580 700 750 750 710 
Hour 13 14 15 16 17 18 19 20 21 22 23 24 
Load(MW) 690 700 730 720 680 610 630 660 620 560 480 430 
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